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1 Introduction 



In Einstein's General Relativity (GR) there are two fundamental physical parameters, the gravitational coupling G 
and the cosmological term A( 4 ), which are typically assumed to be constants. However, the data emerging from the 
experimental and/or observational verification of this assumption arc not conclusive. As a matter of fact, recent 
measurements point toward the possibility that small variations of these parameters cannot be excluded a priori ^ . 

The accuracy in the determination of G is rather poor. Indeed, in the best case it does not exceed 10 -4 . Besides, 
some experiments give inconsistent values for G, which do not overlap within their range of accuracy. Taking into 
consideration the progress and increase of precision in measuring technology, this can be a sign of possible range 
variations of G which might be induced by some "new physics" and/or non-Newtonian interactions [2], [3]. 

In astronomical and geophysical experiments, measurements of possible variations of G with cosmic time provide 
experimental bounds on G/G, which are not very tight. They span from 10~ 10 to 10~ 12 yr _1 depending on the 
experiment and/or observation. Therefore the time variation of G in terms of the Hubble parameter H can be 
written as 

§ = (1) 

where g is a coefficient whose observational bound is |<?| < 0.1. A comprehensive and updated discussion of the 
various experimental and observational constraints on the value of g (as well as on the variation of other fundamental 
"constants" of nature) has recently been provided by Uzan 4 . Therefore, here we will not review the cosmological 
constraints on g imposed by the cosmic microwave background and nucleosynthesis. Instead, we will refer the 
interested reader to 0] and references therein. On the other hand, there are a number of theories, such as scalar- 
tensor theories of gravity [S] and multidimensional cosmological models 6 , which lead to estimates for the variation 
of the gravitational coupling similar to QJ. 

The introduction of a cosmological term in GR is usually considered artificial. Therefore, commonly it is set equal 
to zero. However, the possibility of a small but non-zero cosmological term is suggested by some recent observational 
results from Type la supernovae in distant galaxies and the age problem 0, |E]- Present data reveal that the energy 
density of the field (or quintessence) associated with A( 4 j exceeds the density of ordinary matter [3]. There is an 
extensive literature suggesting that the relation 

A (4) ~ff 2 , (2) 

plays a fundamental role in cosmology HJ. This relation has been obtained in a number of empirical models |1U| . 
|11| . More recently, it was obtained from a model based on the quantum gravitational uncertainty principle and 
the discrete structure of spacetime at Plank length |13| - The dependence Am ~ H 2 explains the current 

observations successfully and provides a much needed large age of the universe. 

In a recent paper |14j we showed that the introduction of a time- varying (744 in brane-world theory yields a number 
of cosmological models which have good physical properties but do not admit constant values for G and Am. In 
particular we considered the five-dimensional metric 

dS 2 = y 2 dt 2 - t 2 ' a y 2 l^- a \dr 2 + r 2 (d9 2 + sin 2 Odtf)] - a 2 (l - a)- 2 t 2 dy 2 , (3) 

where a is a constant. This is a solution of the Einstein field equations in five-dimensions |15| . In four-dimensions (on 
the hypersurfaces y = const.) it corresponds to the 4D Friedmann- Robertson- Walker models with flat 3D sections 1 . 
We used © as the generating 5D space for a Z 2 symmetric brane-universe filled with perfect fluid, with energy 
density p and pressure p satisfying the equation of state p — jp. From Israel's boundary conditions we found that 
the effective G and A(4) vary as 

G fcfm [3(7+1) -«] 

aH, 8irG= (5) -H, (4) 



',2 

G 3(7 + 1) 



and 

_ [3(7+1) -a] 2 2 
A (+ " J^TW ' (5) 

1 The energy density p e f j and pressure p e f t of the effective 4D matter satisfy the equation of state p e j t = np e tf, where n = (2o/3 — 1). 
Thus for a = 2 we recover radiation, for a = 3/2 we recover dust, etc. 
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respectively. Thus, the metric (0 leads an exact brane-model for a FRW universe where where (QJ and J2J hold. 

A crucial point here is that in brane-world theories the quantities G and A( 4 ) are related to each other through 
the vacuum energy density of the 3-brane. They are not independent, as in Jordan-Brans-Dicke and other multidi- 
mensional theories. Thus, either both are truly constants or they vary simultaneously, viz., 

A (4) ~ GG. (6) 

This places particular constrains upon the evolution of brane-world models. Namely, when we take into account the 
experimental bounds (JIJ and we obtain an equation for the Hubble parameter which has to be solved together 
with the generalized Fricdmann equation on the branc. 

Another important feature here is that a cosmological term in 51? induces "natural" constant scales in time, say 
t s , and length (l s = ct s ) in the physics on the brane. Therefore, the evolution of the brane universe can be separated 
into three epochs. For t << r s , the evolution of H and A( 4 j does not differ much from the familiar H ~ t^ 1 and 
Am ~ t~ 2 , in agreement with the natural dimensions of these quantities in AD. For t ~ r s , the overall evolution 
of the expansion scale factor of the brane is dictated by the constant scale. For t >> r a , the evolution depends on 
whether the universe is embedded in a de Sitter or anti de Sitter five-dimensional bulk. 

In this work we investigate the consequences of Q and J5J in the framework of brane-world models. The natural 
question to ask is whether © could lead to a generic analysis of the time variations of A( 4 j and G. Unfortunately, 
the answer to this question seems to be negative. We will see that in order to be able to integrate the equations we 
need to introduce some additional hypothesis or assumption. 

Regarding G, here we consider the simplifying assumption that g is a non-vanishing constant (otherwise G and 
A(4) are constants). Beforehand, we should note that this assumption is not equivalent to requiring G ~ H, in 
general. The physical meaning of this assumption is that the variation of g is much "slower" than that of H and G, 
namely, \g/g\ << \H/H\, \g/g\ << \G/G\. Besides, the numerical value of g should be small enough as to ensure 
that one obtains the correct abundances, in accordance with the discussion in g], and not contradict nucleosynthesis. 

Regarding A( 4 ), without any loss of generality we can write 

A(4)=£W# 2 , (7) 

where £ is some function of time. This form of A( 4 ) is convenient for our purposes and is consistent with the fact that 
supernovae data give us A( 4 ) ~ H 2 at present. In general £ is not constant, but a function which can be obtained 
from the equations of the model. We will see that, in models with non-vanishing bulk cosmological constant, A( 4 ) is 
not strictly proportional to H 2 , i.e., A( 4 ) ^ const x H 2 . These two quantities become exactly proportional to each 
other only asymptotically in time. This is an interesting feature specially in view of latest supernovae results 
suggesting that A( 4 ) was less important is the past and therefore not proportional to H 2 . 

The purpose of this work is twofold. Firstly, to determine the general evolution of the expansion scale factor, 
as well as, the variation of G, A( 4 ) and H in a brane-universe which is compatible with and @. Secondly, to 
examine in some detail the effects of a five-dimensional cosmological constant on the evolution of the brane universe, 
and perhaps on the signature of the extra dimension. 

Our analysis is rather universal since we make no reference to any particular solution of the field equations in 
5D. For generality, we do not impose the signature of the extra dimension either. We consider different possible 
scenarios where the bulk cosmological constant can be zero, positive or negative. For these scenarios, we show that 
models with spacelike extra dimension agree with the observed accelerating universe. 

The extra dimension can be timelike only if the cosmological constant in the bulk is positive. The corresponding 
cosmological models are well behaved and exhibit interesting physical properties, but they do not seem to be of much 
observational significance. 

This paper is organized as follows. In Section 2 we give a brief summary of the theory and of the generalized 
Friedmann equation on a spatially homogeneous and isotropic brane. In Section 3 we show how to incorporate a 
varying vacuum energy into the scheme. We analyze the compatibility of the generalized Friedmann equation with 
the observational requirements and (J2J) • We find that the coupled evolution of G and A( 4 ) completely determines 
the behavior of the brane, except for some adjustable parameters which can be evaluated using some observational 
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data. In Sections 4, 5 and 6 we present a detailed study of the behavior of the brane-universes under consideration. 
In Section 7 we present a summary and discussion 



2 Field equations in bD 

In order to facilitate the presentation, and set the notation, we give a brief review of the equations in brane-world 
theory. We consider the metric 

dS 2 = g liv {x p ,y)dx> 1 dx v + e&{x p ,y)dy 2 , (8) 

where e = — 1 or e = +1 depending on whether the extra dimension is spacelike or timelike, respectively. 
Everywhere we use signature (H e). The Einstein equations in five dimensions are 



^G AB = ^R AB - \g AB ^R = k 2 { ^T AB , 



(9) 



where kr 5 ) is a constant introduced for dimensional considerations and ^T AB is the five-dimensional energy- 
momentum tensor. 

These equations contain the first and second derivatives of the metric with respect to the extra coordinate. These 
can be expressed in terms of geometrical tensors in 4D. 

In absence of off-diagonal terms (g^ = 0) the dimensional reduction of the five-dimensional equations is par- 
ticularly simple |17| . The usual assumption is that our spacetime is orthogonal to the extra dimension. Thus we 
introduce the normal unit (n A n A = e) vector, orthogonal to hypersurfaces y = constant, 



n A = ^, n A = (0,0,0,0,e$). 



Then, the first partial derivatives can be written in terms of the extrinsic curvature 



(10) 



(11) 



The second derivatives, (d 2 g fJlV /dy 2 ), can be expressed in terms of the projection ^C^i of the bulk Weyl tensor 
in five-dimensions, viz., 



(5) 



C A 



BCD 



- (5) 



2 1 

Rabcd - ^ 5) Ra[c9d]b ~ {5) Rb\c9d}a) + g (5) R9a[c9d\b- 



(12) 



The field equations © can be split up into three parts. In terms of the above quantities, the effective field equations 
in AD are, 



(4) G 



h 2 

3 fc (5) 



{5) T aP + {^Tt~\^T)g aP 



e (K a \Kp - KlK aP ) + -g aP (K Xp K xp - {K x ) 2 ) - eE af3 , 



(13) 



where 



E a p = {5) C aAl3B n A n B 
1 dK a0 



+ K ap Kt 



a;/3 



k 2 
. ft (5) 



(5) T aP + (< 5 >T 



\- l f ] T)g aP 



$ d y • ~ ap "P $ 3 
Since E pv is traceless, the requirement E p = gives the inhomogeneous wave equation for <E>, viz., 



(14) 



dK 



= -e— - <S>(eK Xp K Xp + ^ R 
dy 



4)- 



(15) 
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which is equivalent to = kf^S^T^ from ©. The remaining four equations are 



"(5) 

^(^-^)=fcf 5) ^. (16) 
In the above expressions, the covariant derivatives are calculated with respect to g a p, i-e., Dg a p = 0. 

2.1 The brane- world paradigm 

In the brane-world scenario our space-time is identified with a singular hypersurface (or 3-brane) embedded in an 
AdS§ bulk |18j . i.e., it is assumed that the five-dimensional energy- momentum tensor has the form 



(5) 



T AB = A (5) 3AB, (17) 



where A( 5 ) < is the cosmological constant in the bulk. For convenience, the coordinate y is chosen such that the 
hypersurface £ : y — coincides with the brane. Thus, the metric is continuous across E, but the extrinsic curvature 
K^ v is discontinuous. Most brane-world models assume a Z 2 symmetry about our brane, namely, 

dS 2 = g^(x p ,+y)dx^dx l/ + e<5> 2 (x p , +y)dy 2 , fory>0 

dS 2 = g^{xP,-y)dx p dx v +€<$> 2 {xP,-y)dy 2 , fory<0. (18) 

Thus 

|s+= |e- • (19) 
Therefore the field equations in the resulting Z2-symmetric brane universe can be written as 

( 5 )^ = ^(A (5) ^ + ( 5 'f;; H ), (20) 

where g~AB is taking as in 

m and (5)T™, with WfZ } n A = 0, is the energy-momentum tensor of the matter 

on the brane 

(5, 4 ttRe) =^^vf. (21) 
The delta function expresses the confinement of matter in the brane, hence 

p /2 

/ ^ I fx\ -(brane) 

t„ u ( X p,0) = lim / ^T" '<$>dy. (22) 

From Israel's boundary conditions 

&p,v |E+ -K-nv |E-= ~ ek (5) \jl*v ~ 

and the Z 2 symmetry 



| E + -Kpv | s -= -efc (5) t m „ - -g^r , (23) 



we obtain 



K^u |e+= Ie-= ~f fc (5) ~ ^.9^ r ) ■ (24 ! 

2e , 

= _ TT~ - • ( 25 ) 



Then from {El and {EJ it follows that 

- 0. (26) 

Thus represents the total, vacuum plus matter, conserved energy-momentum tensor on the brane. It is usually 
separated in two parts, 

T llv = (Tg llv + T llv , (27) 
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where a is the tension of the brane in 5D, which is interpreted as the vacuum energy of the brane world, and T^ v 
represents the energy-momentum tensor of ordinary matter in AD. From (|24(l . I|25|l and (|27|l we finally get 

<*(5) (- 1 



K„u | E+ = f- \ Tfj.v - pU T + ff )J • (28) 

Substituting Ij28(l and (|17J) into Ijl3|) . we obtain the Einstein equations with an effective energy-momentum tensor in 
AD as 

( 4) G^ = A (4) .g^ + 87rGT^-e^ 5) n^-e£;^, (29) 

where 

A(4) = ^ 5 )( A (5)- £ ^)> (30) 
8nG = -e^-, (31) 

and 2 

= - ^TT^ - \g» v T a0 T<*P + ^g^T 2 . (32) 

All these four-dimensional quantities have to be evaluated in the limit y — > + . The above constitute the basic 
equations in brane- world models. They contain higher-dimensional modifications to general relativity. Namely, local 
quadratic energy-momentum corrections via the tensor \A^ V , and the nonlocal effects from the free gravitational field 
in the bulk, transmitted by E^ v . Another important novel feature is that they provide a working definition of the 
fundamental quantities A(4) and G. 

2.2 Cosmological settings 

In cosmological applications the five-dimensional metric (JSJ) is commonly taken in the form 

dr 2 



dS 2 = n 2 (t,y)dr - a 2 (t,y) 



(1 - kr 2 ) 



+ r 2 (d6 2 +sin 2 



e<5> 2 {t,y)dy 2 , (33) 



where k = 0, +1,— 1 and i,r, 9 and </> are the usual coordinates for a spacetime with spherically symmetric spatial 
sections. 

The metric coefficients are subjected to the conditions 

n{t,y)\ brane = l, a(t,y)\ brane = a (t). (34) 

In this way the usual FLRW line element is recovered on the brane with a$ as scale factor. 

The corresponding field equations in the bulk can be written in a very compact form in terms of the function F, 
which is a first integral of the field equations |24) . namely, 

. . , , (da) 2 (a'a) 2 . „. 

F{t,y) = ka 2 + K -J r + e^, (35) 

where a prime denotes a derivative with respect to y. The (§) and (\) components of © become 

F' = ^ 5) (5) T°, (36) 

and 

F = ^ 5) (5) Tt (37) 

2 With this choice of signs, for perfect fluid IT^,, = (1/I2)[p 2 u fl u l , + p(p + 2p)/i M „] where /i M „ = u^Uv — guu- 
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Now, the field equations (\) = (J) = (|) reduce to 

_ *f 5 ) A (5) ( a'\ 2 k C 



a 



(38) 



s na J 6 V a< ^7 fl2 fl4 

where C is a constant of integration. 

The ordinary matter on the brane is usually assumed to be a perfect fluid 

T l _ lv ^{p + p)u^u lJ -pg^ u , (39) 

where the energy density p and pressure p satisfy the isothermal equation of state, viz., 

P = 7P, < 7 < 1. (40) 

Thus, the boundary conditions J23J), the Z 2 symmetry if24"fl. and (J2SJI yield [HI 



2r 



P(t) "^ } (7 -D<T> 



| brajie 

and 



(41) 

brane 



2c 



. a n 
3 7 + 2 - + - 
a ft 



(42) 

6rane 



These equations show that the tension of the brane and the energy density depend on the details of the model. They 
enable us to investigate the very intriguing possibility that a, and consequently G and Am, might vary with time. 
Evaluating l|38|) at the brane we obtain the generalized Friedmann equation 

^) 2 =A (4)+ 8. G p-^-§ + 5. (43) 

Except for the condition that n = 1 at the brane, this equation is valid for arbitrary $(t, y) and n(t,y) in the 
bulk. This equation allows us to examine the evolution of the brane without using any particular solution of the 
five-dimensional field equations. In what follows we will omit the subscript 0. 

3 The equations for varying vacuum energy 

In the above equation G and A(4) are usually assumed to be "truly" constants. However, as we have already discussed 
there are several models, with reasonable physical properties, for which a variable $ induces a variation in the vacuum 
energy a In this section we show how l|43l) should be modified as to incorporate, or accommodate, the variation 
of these fundamental physical "constants" , matching observational predictions. 
From J2HI) and it follows that 

For perfect fluid l|39l) this is equivalent to 



*,v + T»=0. (44) 



p+(p + p)<d = -&, 
{p+p)a v +p,\hl = a,\hl, (45) 

where = is the expansion, a„ = Uy t \u x is the acceleration and h^ v — u^u v — g^ u is the projector onto the 
spatial surfaces orthogonal to u^. In homogeneous cosmological models the second equation above is empty; only 
the first one is relevant. Namely, for the case of constant vacuum energy a, and the equation of state p = jp, it 
yields the familiar relation between the matter energy density and the expansion factor a, viz., 

(46) 

For the case where the vacuum energy is not constant we need some additional assumption. For example, if a is a given 
function a = a (a), then we integrate the conservation equation (|45l) and substitute the resulting function p = p(a) into 
(|43|l . and thus obtain the corresponding Friedmann equation. However, we have so far no theoretical/observational 
arguments for the evolution of a in time. 
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3.1 G/G-H 

The time variation of G is usually written as (G/G) = gH , where g is a dimensionless parameter. As we have already 
mentioned in the Introduction, nucleosynthesis and the abundance of various elements are used to put constraints 
on g. The present observational upper bound is \g\ < 0.1 In what follows we assume that g is constant. 

Since G ~ a and H = a/a we have a (a) = foa 9 , where fo is a constant of integration. Thus, 

p + 3p( 7 + l)- = -foga^a. (47) 
a 

First, we consider the case where p can be expressed in a way similar to Q46[). i.e., as a power function of a. Therefore, 
we find g = -3D (7 + l)/(/o + D) and 

p = Da- 3 ^ +1 >/^ +D \ (48) 
where D is a positive constant. In order to simplify the notation we set fo = FqD and 

which is ft — (7 — Fo)/(l + F ) or F = (7 — + !)• With this notation we have 

- D - D ^~P) (w\ 

P 3(/S+l)' CT (^ + l) a 3(/3+D- ^ 

Notice that Fo 7^ ((3 ^ 7), otherwise G = 0. We also have, 

| = -3(/?+l)iJ. (51) 

Before going on, we should become aware of the observational bounds of (3. The lower bound comes from the obvious 
requirement dp/ da < 0, while the upper one comes from the observation that \g\ < 0.1 Thus, 

-1<P< -0.966, 0= -0.983 ±0.016. (52) 

We will see that (3 is related to q and Cl p , the deceleration and density parameters, respectively 3 . 
In order to obtain the evolution equation for a we substitute ()50|l into (|43|l 



Mr) -— a (5)-^-^tjttJ ^ny-^ + ^i- ( 53 ) 



Let us now introduce the quantities 



a 3(/3+D 



A = -l^l( 7 + i) 2 ^, 



3fc 2 

C = ^(/3 + l) 2 A (5) , (54) 



in terms of which (|53(l becomes 



(^j 2 = A + Cx 2 - 9 09 + l) 2 [fcx( 6 ' 3+4 )/( 3 ' 3 + 3 ) - Ca; ( 6/3+2 '/( 3/3 + 3 '] . (55) 



^Because of the difficulty of reliably determining g, the limit set by 1521 is just exploratory. We use this number here, not to make 
specific observational suggestions, but in order to be able to discriminate between what is possible or not. 
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This equation admits exact integration, in terms of elementary functions, for a wide variety of parameters A, C, 
13 and C. However, we are interested in the physical models for which A4 ~ H 2 . This introduces some additional 
features which we discuss bellow. 

Let us now notice that l|48|) is not the general solution of ijUJl, which is 



V 



P = 



kg 



a 3(7+i) (3 7 + 3 + g) 



(56) 



where T> is the integration constant. Indeed, setting T> = and denoting g = — 3(/3 + 1) and D — /o(/3 + l)/(7 — (3), 
we recover the solution (|50J) . We will see that, for the models under consideration, the observational requirement 
demands V = 0. 



3.2 A (4 ) ~ H 2 

Without loss of generality we now set 

A( 4 ) =£tf 2 , 

where £ is an unknown function of time which we will determine bellow. Then using (|30|) and (| 5 f )|> we get 



1 2 

Am = -fc (5 )A(5) 



fc 4 5)J D 2 ( 7 -/3) s 



H4) g ^12(/3+l)2 a 6(/3+l) 

In terms of the notation 154|l , this equation becomes 



/dx s 









A 



(7 + I) 2 



■Cx* 



Thus, equating (|53|) and l(5^|) we get 

"3(7 -P) 2 



A 



£(7+ I) 2 



^(6/3+4)7(3/3+3) _ Ca ,(6/3+2)/(3/3+3) 



0. 



(57) 



(58) 



(59) 



(60) 



This equation should take place for all t and x. Therefore, it is important to notice that for the allowed values of /?, 
given by l|52|) . the powers of x in front of k and C are far from being zeroth or second order. Consequently, from (I6U|) 
we get two compatibility conditions. The first one is a condition of physical nature. Namely, 



k = 0, C = 0. 



(61) 



Thus, the resulting cosmological models have flat space sections. This is compatible with astrophysical data 4 from 
BOOMERANG gS] and WMAP Besides, the constant C, which is related to the bulk Weyl tensor and 

corresponds to an effective radiation term, is constrained to be small enough at the time of nucleosynthesis and it 
should be negligible today |27|. 

The second condition consists of three different options of mathematical nature: 



1. If A ^ 0, C ^ 0, then © requires f = 3 and (7 - f3) 2 

2. If C = 0, then 



(7 + 1) 2 ; 

3(7 -P? 



(7+ 1) 2 



(62) 



3. £ is a function of £. 



4 In order to avoid misunderstanding, here do not use cosmological parameters emerging from BOOMERANG or WMAP. The point is 
that the conditions in 1611 . for a model which does not obey the usual Einstein equations, do not contradict well known results obtained 
under the assumption that the (usual) Einstein equations are valid. 
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Certainly, the first option cannot take place by virtue of Q52|): only the second and the third are possible ones. We 
will see them at work in Sections 4 and 5, respectively. 

If we apply the same arguments used above to the general solution (|56|l and a = foa 9 , instead of i|50|) . we get two 
compatibility conditions. The first one for T> ^ is very restrictive and requires g = — 3(7 + 1), which we disregard 
on observational grounds since \g\ < 0.1. The second condition for T> = yields the case discussed above. 

The conclusion from this section is that the cosmological models where G/G ~ H and A( 4 ) ~ H 2 have flat space 
sections and zero (or negligible) Weyl radiation from the bulk. 
Their time evolution follows from 

/ot- 11 -* <63) 

where to is a constant of integration. 

The solution of this integral depends on the constants A and C, which contain e and A/ 5 ). Thus, it allows us to 
study some important questions such as: How does the bulk cosmological constant affect the behavior of the brane 
universe? Does it influence the choice of the signature of the extra dimension? 

In the next Sections we analyze the solutions of l(B3*j l for Ar 5 \ = 0, A( 5 ) > and A( 5 ) < 0. 



4 Bulk with A (5) = 

We consider first the behavior of the model with vanishing bulk cosmological constant. This case is important 
because, as we will see later, it is the limit of the models with A( 5 ) ^ and spacelike signature. 

From l|54[l it follows that the extra dimension has to be spacelike (e = —1). The expansion factor is given by 



o(t) 



1/3(0+1) 



(64) 



where we have set to = 0, in such a way that the big bang occurs at t = 0. Here (3 is related to the deceleration 
parameter q — —da/a 2 as 

q = 2 + 3(3. (65) 

Therefore, the ratio 1)51(1 becomes GjG = —(q+ l)H, which is the same as in separable models jUj. For the vacuum 
energy we find 

2 (7 - P) 



Consequently, the 4D effective cosmological term Au) and the gravitational coupling G vary as 



A <o= ( 7 + i) 2 



where H = a/a is the Hubble parameter, viz., 



H 



8irG 



1 1 



^ 5 ) (7-/3) 



7+1 



-H, 



3(0+1) t' 



(66) 



(67) 



(68) 



Thus, here the second option (|62|l is taking place. There is a remarkable connection between (3, 7 and the density 
parameter fl p . Specifically, 



^G P _ _ 2( 7 -/?)(/3 + l) 
3H 2 p (7 + I) 2 



(69) 



Thus, giving 7 and Q p we obtain f3, viz., 



(/3+i) = -(7+i)(i±vr 3 2o;j ) 



(70) 
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which gives an upper bound for Q p . Namely, 

< \ (7i) 

regardless of the specific value of 7 and q. 

From the above equations, it follows that the density parameter Q p can, in principle, be determined by measure- 
ments of the deceleration parameter q 

_ 2(3 7 -g + 2)( g + l) 

This is a useful feature which results from the assumption that A/ 5 ) vanishes. A similar situation occurs in the 
familiar FRW models. We also note the relationship between the density parameters, viz., 

a (4) _ ggq + TO; 

~2~jj2 ^ 4(1 + g)a • ( 73 ) 
4.1 Behavior of the model for different Vt p 

Let us study the behavior of the model for different values of the density parameter fi p , within the range of values 
allowed in (J7TJ • From and (J7UJ) it follows that 



a(t) ^ t 2/3( 7 +i)[i±yT32?y (?4) 

Thus, in the upper limit, for fi p — > 1/2, the evolution of the scale factor of the universe matches the one in familiar 
FRW models, viz., a(t) « i 2 / 3 (T +1 ). Accordingly, the deceleration parameter (|65[1 reduces to q w (1 + 37)/2 which is 
the usual one in FRW cosmologies. An interesting feature of the model is that, in this limit neither the gravitational 
coupling is constant, nor the cosmological term becomes zero. Instead we have 

^G=^H, A (4) = ^, H=^\, (75) 

where H is the usual Hubble parameter in FRW models, as expected. 
Let us now study the behavior of our model for other values of Q p . 

It is clear that the deviation from the FRW models increases as fi p moves away from 1/2. But for every value of 
f2 p there are two possible models; one for each sign in front of the root in l|7()|l. They both reach the FRW model for 
fip = 1/2. The question now is whether they satisfy physical conditions. 

Since a reliable and definitive determination of Q p has thus far eluded cosmologists, in our discussion we consider 
several values of Q p , although fl p w 0.1 — 0.3 seem to be the most probably options. 



4.1.1 Accelerated expansion 

We consider the behavior of dust-filled universes (7 = 0). In Table 1 we illustrate the evolution of such universes 
corresponding to the negative sign in (|70|l . It is interesting to note that there are no arbitrary parameters or constants 
in the solution. Specifying fi p we find (3, q and Qa = A( 4 )/3i/ 2 . The age of the universe is calculated assuming 
JJkOJxIQ-V 1 - 



Table 1: A( 5 j = 0. Accelerated expansion 




n A 




<?(-) 


(^_)/10 lu ) yr 


0.5 


0.25 


-0.500 


0.500 


0.95 


0.4 


0.522 


-0.723 


-0.170 


1.71 


0.3 


0.664 


-0.816 


-0.448 


2.58 


0.2 


0.783 


-0.887 


-0.661 


4.21 


0.1 


0.889 


-0.947 


-0.841 


8.98 


0.08 


0.907 


-0.958 


-0.874 


11.33 


0.04 


0.936 


-0.979 


-0.938 


22.67 
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This case is interesting because it shows an accelerated expansion of the universe, in agreement with modern 
observations. The acceleration is driven by the repulsive effect of the "dark energy" associated with A( 4 ) , which 
clearly dominates the evolution here. We notice that f2 p + ^ 1. There is a contribution, fl p 2, from the quadratic 
correction in the generalized FRW equation l|43|). so that £l p + SI a + f2„2 = 1. This contribution decreases as the 
universe gets older. For £l p < 0.4, the universe is quite old and tt p 2 is "negligible". 

We have already mentioned that for fl p — 0.5, the expansion scale factor is the same as in the de Einstein-de 
Sitter solution, although the matter content is totally different. We now see that Vt p 2 = Sl^ = 0.25 for this solution. 
This is another example of the well known fact that in GR the same geometry can be attributed to different matter 
distributions. 

Dynamical mass measurements from WMAP Mission reveal that the matter content of the universe is about 27% 
of the critical density (4% Atoms, 23% Cold Dark Matter). The rest 73% is usually declared either as "missing", or 
Dark energy. Also, according to recent measurements 5 the acceleration parameter is, roughly, —0.5 ± 0.2. We note 
that the entries in the third and fourth rows are consistent with this picture of the universe. Another remarkable 
feature here is that the age of the universe is much larger than in the usual FRW model (first row). 

4.1.2 Decelerated expansion 

We now consider the solution with positive sign in front of the root in (|70(l . For the purpose of comparison, we again 
consider the evolution of dust universes under the same set of values of Q p as in the model with negative sign. The 
relevant parameters are presented in Table 2. 



Table 2: A( 5 ) = 0. Decelerated expansion 


n p 






<?(+) 


(T i+) /10^) yr 


0.5 


0.25 


-0.500 


0.500 


0.95 


0.4 


0.076 


-0.276 


1.170 


0.65 


0.3 


3.37 x 10"^ 


-0.184 


1.448 


0.58 


0.2 


1.27 x 10~ 2 


-0.113 


1.661 


0.53 


0.1 


0.37 x lO"^ 


-0.053 


1.841 


0.50 


0.08 


0.17 x 10~ 2 


-0.042 


1.874 


0.49 


0.04 


0.41 x 10" 3 


-0.020 


1.938 


0.48 



In the present case the repulsion associated with the cosmological term is negligible. Here the evolution is 
dominated by the quadratic correction term. As an illustration consider the entries on the third row, for which 
f2 p 2 ss 0.7. The huge gravitational attraction produced by this term, for a spacelike extra dimension, explains the 
large deceleration parameter. 

This solution presents a number of interesting features. For example, a decrease in f2 p entails a decrease in the 
age of the universe. This is the opposite of what we see in Table 1. Also, T( + ) is almost the same for H p w 0.1 — 0.3, 
while in Table 1 the age of the universe changes almost four times in the same range. 

However, the solution with positive sign in l|70l) seems to have little in common with present observations. Firstly, 
the universe is quite young. 6 Secondly, it does not fit the observational requirements on the ratio G/G, i.e. on (3, 
neither on the deceleration parameter. 

We would like to finish this Section by stressing the fact that the parameter f3 is related to (i) the ratio G/G, (ii) 
the deceleration parameter, and (iii) the density parameter Q, p through the equations Q52JI . I|65(l and (|70|1 . respectively. 

5 We would like to emphasize that here we are only using the measured values of q, the rest of the parameters are obtained from the 
model. 

6 This is way the quadratic corrections £l p 2 dominate. 
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5 de Sitter bulk, A( 5 ) > 

Although in brane-world theory our universe is embedded in a higher-dimensional space with negative cosmological 
constant, the solutions to the evolution equation (|63[) depend analytically on A( 5 ), allowing us to continue to the 
range of positive bulk cosmological constant. In this Section we explore the case with A( 5 ) > 0. An attractive feature 
of this case is that the extra dimension can be either spacelike or timelike. We will discuss these two cases separately. 

5.1 Spacelike extra dimension e = — 1 

The integration of l|63[l yields 

(/3 + l)V 6A (5) 

where once more we have set to — and 



C = %)(/<' - "\/^- I") 
The deceleration parameter here is 

q = 2 + 3/3- 3(/3+ l)tanh 2 \/Ci. (78) 

The matter and vacuum density vary as 

2VC 2(7-/?) VC 

o= =-, u — 79 

kf 5) (7 + l)siuh y/Ct k 2 {5) + ([3 +l)sinh VCt 

As in the previous solution we have G/G = — 3(f3 + 1)H and 



( 7 + l) 2 -( 7 -/?) 2 



_ 3(7-/3)^ 2 fc (5) A (5) 



(80) 



8.r k MtZ® m , k ^ (R1) 
8 * G - (7 + 1) ~w^' (8 ) 

where 

H= ^ = . (82) 

3(/3+ l)tanh7Ct 

We note that A/4) is always positive and H is a monotonically decreasing function of time bounded bellow by 
H = fc( 5 ) ^/A( 5 )/6. This assures the positivity of G, as long as 7 > 0. 
The general behavior of the solution is as follows 

1. For A(5) = we recover the previous solution. 

2. For small values of t, near the big bang, this model behaves exactly as the previous one Q64 fl -H68 fl . 

3. At large times the universe is expanding with positive acceleration, q w — 1. The expansion is exponential and 
A( 4 ) tends to A = fc 2 5 ^A( 5 )/2. 

Now, if we combine (|78|) and l|82|l we obtain 

C = 3H 2 (2 + 3(3- q) (1 + (3) . (83) 



13 



Using this expression into 8wGp from (|79|l we get 



3ft P (7 + l) 



2(1 + 9) 

The condition (f3 + 1) > sets an upper limit on the density parameter, viz., 

P 3(1 + 7)' 

We also obtain an interesting expressions for Am . Namely, 

"(7-/3) 2 (l + g) + (l+7) 2 (2 + 3/3- <Z ) 



A (4) = 



(l + /3)(l+ 7 ) 2 



if 2 



(84) 



(85) 



(86) 



Coming back to 1)57(1 . here the factor (£) in front of iJ 2 is noi constant, because q is a function of time. Thus, this 
solution corresponds to the third option mentioned after (|62|l . 
For the gravitational "constant" we get 



V3( 7 + i)%/TT7? 

The above equations Q84[l. (|86|l relate the observational quantities 7, f2, g and /3. The age of the universe is 



T 



C 



tanh 



2 + 3/3 - <, 
3(/3 + l) 



1/2 



(87) 



(88) 



where C is obtained from (|83|) . Finally, the five-dimensional quantities /c( 5 ) and A( 5 ) can be evaluated from l|87|l and 
(|77|l . We note that there are no free parameters left in the solution. 

5.1.1 Characteristic time 

It is important to notice that a non-vanishing cosmological constant in the bulk induces a natural time scale, in AD. 

c? 5 )A(5)/2. Thus, t s — 3(f3 + l)/VC or, in terms of observational quantities 

l 1/2 

(89) 



We define it as r s 



1 



H 



4(1 + g) 2 - 9^(1 + 7) 
4(1 + q) 2 n A - 9fi 2 (l+ 7) 2 



We remark that although the quantities on the r.h.s. are functions of time, r s is a "universal" constant fixed by 
the five-dimensional embedding bulk. This equation shows how we can evaluate this constant from measurements 
performed in AD. The denominator vanishes and t s — > 00 when the density parameters are related as in l|73|l. Thus, 
the "upper" bound for t s takes place in the models discussed in the previous Section. 
In order to get an expression for the lower bound, we rewrite l|89|) as 



4(1 + q) 2 _ H 2 T 2 S - 1 

9^2(1 + 7 )2 - n A H 2 T 2 - 1 ' 



(90) 



From l)84|) it follows that the numerator in (|89|l is always positive. Then, from J7a < 1 we obtain (H 2 t 2 — 1) > 0. 
Consequently from lO we get (fl A H 2 T 2 - 1) > 0. Thus 



1 



< T, < OO 



(91) 



For a universe with w 0.7 and H w 0.7 x 10 10 yr 1 , the lower bound for the characteristic scale is rs 1.7 x 10 10 
years or 17 billion years, which is more that the age of the universe according to the data from WMAP Mission. 
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5.1.2 Observational constrains 



Let us now use the observational bound (|52|) . From (|84|l we obtain 

30,(1 +7) 30,(1 + 7 ) 2 

2 1 < q < 2(7 + 0.966) A ' 

where we have assumed (7 + 0.966) > 0. This equation implies, q < 0, an accelerated expansion of the universe in the 
range of values allowed for O p . However, not all q and O p generate adequate physical models. We have to take into 
account that here A( 5 ) > 0. Therefore, from (|54|l it follows that C > 0. Consequently, (|83|1 requires (2 + 3/3 — q) > 0. 
If we apply the observational bounds JS2J), as well as (|S5|l . we get 



Kg<-0.88, Q p <-^-, (93) 
(1 + 7) 



which are even more restrictive than l|92|) . According to the present estimates for q and O p we are nowhere near 
these values. Thus, although the present model is attractive from a theoretical point of view, it seems to have limited 
practical application. 

5.2 Timelike extra dimension e = +1 

For completeness we now consider the case where the fifth dimension is timelike, although these models are contro- 
versial in many respects, including the problems of causality and quantization. We are not going to discuss these 
problems here, instead we refer the interested reader to |28| and references therein. 
In this case the scale factor is given by 



a 



3(0+1) _ fc(5)(7 + l)£> 



03 + 1)76A^ 



cosh VCt, (94) 



where once more we have set to = 0. The expressions for the physical quantities p, a, G, Am and H are formally 
obtained from the solution discussed in Section 5. 1 by changing sinh VCt — > cosh VCt and tanh VCt — > 1/ tanh VCt. 
However, there is a qualitative difference between the solutions. 

1. This is a "bouncing" solution where the universe never collapses to a singularity. The universe is always 
expanding with positive acceleration. 



2. The Hubble parameter is an increasing function of time and is bounded above by fc( 5 )-^/A( 5 )/6. 

3. The vacuum energy must be negative in order to ensure the positiveness of G. Namely, 



G= ^- — = , /3> 7 (95) 

3( 7 + 1)03+1) cosh VCt 



We conclude that the model with a timelike extra dimension seems to be of no observational significance because 
for all t we have q < — 1, which is the opposite to what we expect (q > —1), based on recent measurements. 

5.3 Asymptotic behavior and classical inflation 

The late time behavior of solutions with A( 5 ) > is similar for both signatures, e = — 1 or e = +1, although they 
drastically differ at the origin of the universe. In both cases the expansion factor becomes identical to the one in the 
de Sitter solution, 

a(t) - exp yj (A/3)i = exp(i/r s ), (96) 
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as t » t s . It is clear that we are far from this asymptotic regime. It should be emphasized that this exponential 
behavior does not arise from a "false- vacuum" equation of state p = —p as in inflation, because 7 7^ —1 in these 
solutions. The reason for this is that the assumption G/G = gH is equivalent to the requirement that p and cr 
form a combined fluid with energy density p — p + a and pressure p = (3p. Then, the observational constraint 
— 1 < [3 < —0.966 implies that the combined fluid behaves nearly like a cosmological constant, which dominates at 
late times and thus producing inflation. 



5.3.1 Vacuum equation of state 

For consistency, we should now show that this model is compatible with classical inflation for the equation of state 
of false-vacuum, 7 = —1. Indeed, if 7 — —1, then from l|54|l it follows that ^4 = 0. Thus, the evolution equation (|63|) 
requires C > 0, which entails A( 5 ) > 0. Integrating l|63[l we get x = expVC(t — to), which in terms of the original 
notation is identical to the de Sitter expression (|96|l . 

We conclude this section with the following comments. 

(i) From 185|) it follows that for a every given value of q, there is a range of possible values for f2 p . Similarly, for a 
given f2 p , there is a range of allowed values for q, which is given by (|92() . This is different from the case with A( 5 j = 
where they are related by Ij72(l . Therefore, for A( 5 ) ^ the sole specification of one of these parameters is not enough 
to determine the characteristics of the model. Here we have to specify both, independently. This additional degree 
of freedom is a consequence of the introduction of a non-vanishing cosmological constant in the bulk. 

(ii) Measurements in AD allow to "predict" the value of the five-dimensional constants and A( 5 ). Namely, if 
we measure (f2 p ,g), then (3 follows from l|84[l . The value of J7a is then obtained from (|86|l . Finally, from (|87|l and 
(18911 we evaluate the 5D constants fc^ and r s (and/or A( 5 )) in terms of H~ x , If we have (fl p , ^a), then (3 and q are 
given by the solutions of the system of equations 184|l and Q86JI , then following the same steps as above we get the 
rest of the parameters. 



6 Anti-de Sitter bulk, A( 5 ) < 

This case is important, it corresponds to the brane world scenario where our universe is identified with a singular 
hypersurface (or a three-brane) embedded in an AdS^ bulk. 

For A( 5 ) < the extra dimension has to be spacelike and the evolution of the scale factor is given by 

a 3( W)= %)(7 + l)D ^ 
(/3 + l)V^f 

where once more we have set to — and 



|C| = fc (5) (/?+l)y^i (98) 
In the present case the deceleration parameter is 

9 = 2 + 3 / 3 + 3(/3+l)tan 2 V](7|i, (99) 

where 

\C\ =3H 2 (q-2-3/3)(l+f3). (100) 



This model is formally obtained from the solution in Section 5.1 by making the change */A(5) — > i-y/|A( 5 )|. Thus 
\fC — > i^J\C\. Consequently, (^smh\/C~tj fi— » sin -\/|C|i. However, they are drastically different from each other. 
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The present solution (A( 5 ) < 0) represents a spatially flat but recollapsing universe. The recollapse time T rt 
given by sin yJ\C\T rec — 0, which in terms of the characteristic time t s defined in Section 5.1.1 becomes 7 



T — 

■t-re.c. — 



6(/?+l)' 

We note that for the physical values T rec > 15r s . The age of the universe can be written as 



T 



tan 



3/3 



3(0+1) 



1/2 



(101) 



(102) 



In Table 3 we use this solution to estimate the parameters q, in dust-filled universes, with various values of 
il p . We also evaluate, in units of H~ l , the characteristic and recollapse time, r s and T rec , as well as the age of the 
universe T . 



Table 3: AdS bulk, A (5) < 


n p 


^(A) 


q 


Ht s 


HT rec 


HT 


0.5 


0.495 


-0.236 + 0.013 


0.278 


7.943 


6.569 


0.4 


0.596 


-0.389 + 0.011 


0.311 


9.044 


7.650 


0.3 


0.697 


-0.542 ±0.008 


0.359 


10.760 


8.396 


0.2 


0.798 


-0.694 + 0.006 


0.487 


13.004 


9.246 


0.1 


0.898 


-0.846 + 0.003 


1.005 


20.254 


10.048 


0.08 


0.919 


-0.877 + 0.002 


1.211 


25.996 


11.417 



For a given VL p the values of q, as obtained from l)92|l . are spread over a "small" range. In other words, if we know 
ilp, we get q with a great accuracy, which is up to 1% — 5%. We should emphasize that this precision comes from 
the experimental bounds on G/G. Then, we use the mean value of q to obtain [3 from l|84f) . which we substitute into 
(|S5j> to get Q A . 

We note that the entries on rows two, three and four are similar to those in Table 1. Here the deceleration 
parameter q is a monotonically increasing function of time and changes its sign for [3 less than w —0.66. Therefore, 
for the values allowed by 1)52(1 . the universe initially expands with acceleration (q < 0), then changes sign becoming 
positive and q — > +oo for t — > T rec . 

This behavior is a consequence of the effective cosmological term A( 4 ) which is positive near the big bang and, as 
t — > T rec , becomes 



A 



(4) 



fc( 5) A (5 ) 



1 - 



(7 + I) 2 



(103) 



which is negative for A( 5 ) < 0. 

Finally, we would like to note some similarities between solutions with a spacelike extra dimension e = — 1. 
Firstly, they have an analogous behavior near the big bang, regardless of the value of the bulk cosmological constant. 
Secondly, the gravitational coupling of geometry to matter, G, is divergent for t — > 0. This last feature is of particular 
interest, it suggests that the gravitational interaction was much stronger in the past than at the present time. 



7 Summary and conclusions 

In summary, we have studied the consequences of the conditions G/G = gH, with \g\ < 0.1, and A( 4 ) ~ H 2 on 
cosmological models based on the brane-world scenario. These two conditions lead to the requirement k = 0, C = 0, 
in l|61|) . The same physics is obtained if we assume G/G — gH and k = 0, C = 0. In this case A/ 4 j ~ H 2 is not an 
extra condition but a consequence of these assumptions. However, if we choose to postulate A( 4 ) ~ H 2 , and k = 0, 

7 In this case we use the magnitude of A( 5 j, namely, A = fc? g . |A(g) \/2. Therefore, for the evaluation of r B in 1891 we have to take the 
magnitude of the quantity in parenthesis, otherwise r s would be a complex number. 
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C = 0, then we do not obtain G/G = gH, except for the case where A( 5 ) = 0. In other words, we recover only the 
models discussed in Section 4. 

The study of these cosmological models is attractive for several reasons. Firstly, conditions QJ and are 
motivated by a number of theoretical and experimental observations. Secondly, brane theory provides a framework 
within which to consider the simultaneous variation of G and Au\ , because they are related to the vacuum energy 
which in principle could vary with time. Thirdly, the theory has no arbitrary parameters. 

Although we introduce two "new" parameters, besides the observed g (or /3), viz., A( 5 ) and kr^, they can be 
evaluated by means of quantities measured in AD only. Specifically for A( 5 ) we have 

(2 + 3/?- g )( 7 -/?),/T+^ / H 3 \ 
A(5) = V3( 7 + l)(l + /3 )3/ a {^g) • ( 104 ) 

As an illustration we consider the entries on the third row in Table 3, for which we get A( 5 ) w —8.8 x 10~ 46 gr cm~ 3 
s -1 . The value of can be obtained from l|S7|) as k^ « 26.23 x 10 10 gr^cm 3 s^ 1 . Thus, the r.h.s. of the field 
equations in © is a very "small" number, viz., fcLA( 5 ) w 2.3 x 10 -34 s~ 2 . In terms of the Plank units of time 

tpiank = 10~ 43 s, we get fc? 5 sA(5)/2 ~ |A(4)| « 10~ 120 tpf ank , as one expected. Although this number is small, it has 
a crucial influence on the global behavior of the universe. 

The model requires the universe to have flat (k — 0) space sections, which is consistent with observations. Also, 
it is interesting that all possible scenarios, except those in Table 2, agree with the observed accelerated expansion of 
the universe and are dominated by the Dark energy associated with Am . 

In addition, the observed value of q = —0.5 ± 0.2 fits the ones predicted by the models with A( 5 ) < 0, as shown 
in Tables 1 and 3 for tt p w 0.1 — 0.3. This is different from the models with positive bulk cosmological constant. 
These present good physical properties but do not appear to be compatible with current observations, regardless of 
the signature of the extra dimension. 

In all cases the cosmological models with spacelike extra dimension have similar behavior near the big-bang. If 
A(5) = 0, then the induced cosmological term in AD is positive and the brane universe expands continuously in a 
power-law time-dependence. If A( 5 ) > 0, then the universe becomes dominated by a positive cosmological term A( 4 ), 
which tends to a constant value. The effect is an asymptotic de Sitter expansion, which occurs without a false- vacuum 
equation of state and regardless of the signature of the extra dimension. If A( 5 ) < 0, then A( 4 ) is initially positive, 
but later changes its sign. The universe becomes dominated by a negative cosmological term that causes the universe 
to recollapse. 

Consequently an inevitable conclusion of our work is that the universe must recollapse at some time in the future 
if it is embedded in an Anti-de Sitter five-dimensional bulk, which is the usual case in brane models. Fortunately, 
according to Table 3, we are nowhere near the time of recollapse. We emphasize that this behavior is independent 
on the (small) size of the (negative) bulk cosmological constant. 

However, we note that the entries in Tables 1 and 3 are not very different from each other for fl p « 0.1 — 0.3. 
This is explained by the "negligible" difference between the numbers in the r.h.s. of the field equations @ in both 
cases. The entries in Table 1 are calculated for k^A^ = 0, while the ones in Table 3 for kf 5 \A^ w 2.3 x 10~ 120 

tp, k . Such a small difference does not significantly alter the evolution of the universe. 

Therefore, the solution with accelerated expansion of Section 4.1.1 gives a good description, to a first order 
approximation, of the evolution of a AD universe embedded in a five-dimensional spacetime with non- vanishing A( 5 ) 
and spacelike extra dimension. We would like to emphasize that the whole analysis in this paper is independent of 
any particular solution used in the five-dimensional bulk. This is a great virtue of brane-world models as noted at 
the end of Section 2. 

However, one could still ask whether the brane metrics described here can be embedded in a five-dimensional bulk. 
The answer to this question is positive. Clearly the model with A( 5 ) = 0, which generalizes the familiar power-law 
solution characteristic of flat FRW universes, is embedded in the cosmological metric with separation of variables © 
which is discussed in Regarding the models with A( 5 ) ^ 0, they can be embedded in five-dimensional "wave- like" 
cosmologies of the type discussed in |2H) - If m equation (38) of |2E] we take variable a as here in H50JI ■ then the scale 
factor a for such wave-like models is governed by an equation which is identical to H53|) in this paper. 
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We would like to finish with the following comments. The issue of constraining brane universes with observational 
data like supernova, CMB and quasars has been done in detail for brane cosmology in |29| and |30| . Our aim in 
this work has been to study the simultaneous variation of G and A( 4 ) and constraining the model with observations. 
Such study has not been done before. Besides we have very little indication of the 

evolution of G and A( 4 ) in time. Our working hypothesis, that G/G = gH, is an extrapolation of the present limit 
of \G/G\. It is not clear whether it remains valid as early as at nucleosynthesis. Many other important questions 
remain open. Among them the behavior of perturbations and structure formation. Also the issues of obtaining the 
4-dimensional Newton law as well as the quantization in models with positive bulk cosmological constant and models 
with timelike extra dimension, respectively. These questions are beyond the scope of the present work. 
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